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SECTION  A 


An  understanding  of  the  frequency  dependence  of  decoupling  is  essential  for  monitoring 
nuclear  explosions  detonated  in  large  underground  cavities.  Analysis  of  available  data  from  the 
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recorded  at  a  common  station,  KN-UT.  The  spectral  ratio  Pn/Lg  varied  strongly  with  scaled 
depth,  mainly  due  to  significantly  greater  dependence  on  scaled  depth  of  the  spectra  of  Pn  than 
of  Lg.  These  results  are  in  agreement  with  those  from  Salmon  and  Sterling  with  scaled  depths 
of  475  and  1 145  m/(kt)1/3,  respectively.  A  possible  reason  for  the  decoupling  at  higher  frequen¬ 
cies  to  be  different  for  P  and  Lg  is  therefore  the  greater  scaled  depth  of  decoupled  shots.  The 
larger  variability  of  Pn  spectra  than  those  of  Lg  can  be  due  to  an  explosion  (especially  those 
that  are  not  overburied)  approximating  a  pressurized  ellipsoid  cavity  with  major  axis  in  the 
vertical  direction  so  that  the  source  function  is  more  variable  for  down-going  waves  comprising 
Pn  than  for  others  making  up  Lg.  The  spectral  ratio  Pn/Lg  may  be  useful  in  discriminating 
between  decoupled  (or  overburied)  and  normal  shots. 


SECTION  B 

The  amplitude  and  spectra  of  teleseismic  and  regional  seismic  P  waves  depend  on  the 
material  properties  in  the  source  region.  Ihe  effects  of  strain  and  frequency  dependent  attenua¬ 
tion  are  therefore  important  for  understanding  the  source  characteristics  of  coupled  versus 
decoupled  explosions.  Near-field  velocity  and  acceleration  data  from  the  Salmon  and  Sterling 
explosions  in  salt  were  analyzed  for  strain  and  frequency  dependent  attenuation.  Attenuation, 
parameterized  as  1/Q,  was  estimated  for  various  source-receiver  distance  ranges  and  frequency 
bandwidths  by  using  two  methods:  (a)  mean  spectral  ratio  slope  and  (b)  average  amplitude 
spectral  ratio.  Results  from  Salmon  indicate  that  (1)  in  the  frequency  range  of  1-25  Hz,  Q  is 
about  5  to  10;  (2)  Q  increases  with  source-receiver  distance,  suggesting  lower  Q  for  larger 
strain  levels;  (3)  Q  in  the  frequency  range  of  25-50  Hz  is  substantially  higher  than  in  the  fre¬ 
quency  range  of  1-25  Hz,  or  Q  increases  with  frequency.  Analysis  of  data  from  Sterling  indi¬ 
cates  that  Q  is  frequency  independent  and  about  one  order  of  magnitude  larger  than  for  Sal¬ 
mon.  Therefore  the  near-source  attenuation  was  strain  and  frequency  dependent  and  hence 
non-linear  for  Salmon  out  to  a  range  of  600  m. 
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nuclear  explosions  detonated  in  large  underground  cavities.  Analysis  of  available  data  from  the 
nuclear  explosions  Salmon  (5.3  kt  in  salt)  and  Sterling  (0.38  kt  shot  in  the  cavity  created  by 
Salmon) -recorded  at  shot-receiver  distances  of  16,  32,  and  27  km,  confirmed  the  earlier 

findings  of  Blandford  and  Woolson  (1979)  indicating:  (l)  reduced  decoupling  of  Sterling  at 
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higher  frequencies  for  both  P  and  "Lg,"  and  (2)  relatively  greater  decoupling  of  Sterling  for  Lg 
than  for  P  at  high  frequencies.  In  an  attempt  to  understand  why  the  frequency  dependence  of 
the  decoupling  ratio  is  different  for  P  and  the  Lg  group,  we  examined  the  spectral  content  of 
Pn  and  Lg  from  7  Pahute  Mesa,  NTS  explosions  covering  a  wide  range  of  scaled  depths  and 
recorded  at  a  common  station,  KN-UT.  The  spectral  ratio  Pn/Lg  varied  strongly  with  scaled 
depth,  mainly  due  to  significantly  greater  dependence  on  scaled  depth  of  the  spectra  of  Pn  than 
of  Lg.  These  results  are  in  agreement  with  those  from  Salmon  and  Sterling  with  scaled  depths 

of  475  and  1 145  m/(kt)1/3,  respectively.  A  possible  reason  for  the  decoupling  at  higher  frequen- 
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cies  to  be  different  for  P  and  Lg  is  therefore  the  greater  scaled  depth  of  decoupled  shots.  The 
larger  variability  of  Pn  spectra  than  those  of  Lg  can  be  due  to  an  explosion  (especially  those 
that  are  not  overburied)  approximating  a  pressurized  ellipsoid  cavity  with  major  axis  in  the 
vertical  direction  so  that  the  source  function  is  more  variable  for  down-going  waves  comprising 
Pn  than  for  others  making  up  Lg.  The  spectral  ratio  Pn/Lg  may  be  useful  in  discriminating 
between  decoupled  (or  overburied)  and  normal  shots. 
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1  Raw  common-receiver  vertical  component  data  recorded  at  station  10S  4 

(Receiver  location  C2000)  for  the  three  explosions  Salmon,  Sterling, 

and  Stercal.  The  beginnings  of  the  P  and  Lg  wavetrains,  indicated  by 
the  letters  P  and  S,  are  separated  by  about  3.3  sec. 

2  Spectra  of  the  observed  (a)  P  (window  1.02  sec  long)  and  (b)  Lg  (win-  6 

dow  2.05  sec  long)  for  Salmon  recorded  at  10S.  The  spectra  of  noise 
preceding  the  onset  of  direct  P  are  also  shown.  The  amplitude  values 

are  in  arbitrary  units. 

3  Spectra  of  the  observed  (a)  P  (window  1.02  sec  long)  and  (b)  Lg  (win-  7 

dow  2.05  long)  for  Sterling  recorded  at  10S.  The  spectra  of  noise 
preceding  the  onset  of  direct  P  are  also  shown.  The  amplitude  values 

are  in  arbitrary  units. 

4  Spectra  of  the  observed  (a)  P  (window  1 .02  sec  long)  and  (b)  Lg  (win-  8 

dow  2.05  sec  long)  for  Stercal  recorded  at  10S.  The  spectra  of  noise 
preceding  the  onset  of  direct  P  are  also  shown.  The  amplitude  values 

are  in  arbitrary  units. 

5  Theoretical  spectra,  based  on  von  Seggem  and  Blandford’s  (1972)  10 

source  function  for  salt,  for  yields  of  5.3,  0.38,  and  0.0027  kt,  respec¬ 
tively.  The  corner  frequencies  (half  amplitude  points)  for  the  three 
spectra  are  at  about  4,  10,  and  50  Hz,  respectively. 

6  Theoretical  spectral  ratios  for  Salmon/Sterling  (tamped),  Sterling  1 1 

(tamped)/Stercal,  and  Salmon/Stercal  over  the  frequency  range  1  to  100 

Hz. 

7  Spectral  ratio  Salmon/Sterling,  corrected  for  noise,  with  logarithmic  12 

frequency  scale,  for  (a)  P  window  1.02  sec  long,  and  (b)  Lg  window 

2.05  sec  long.  The  dashed  lines  show  mean  (least  squares)  slopes  over 
the  frequency  range  5  to  25  Hz.  The  mean  value  of  slope  and  its  stan¬ 
dard  deviation  (S.D.)  are  indicated. 
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8  Spectral  ratio  P/Lg,  corrected  for  noise,  for  (a)  Salmon,  (b)  Sterling  14 

and  (c)  Stercal  recorded  at  the  station  10S.  Points  for  which  S/N 
power  ratio  is  less  than  2  are  not  plotted.  The  dashed  lines  show  mean 
(least  squares)  slopes  over  the  frequency  range  5  to  25  Hz.  The  mean 
value  of  slope/Hz  and  its  standard  deviation  (S.D.)  are  indicated  for 
frequency  ranges  of  5  to  25  Hz  and  5  to  40  Hz. 

9  Raw  common-receiver,  vertical  component  data  at  station  20S  for  Sal-  16 

mon  and  Sterling.  The  beginnings  of  the  P  and  Lg  wavetrains,  denoted 

by  the  letters  P  and  S,  are  separated  by  about  5.5  sec. 

10  Spectra  of  the  observed  P  arrivals  (window  2.05  sec  long)  for  (a)  Sal-  17 

mon  and  (d)  Sterling  recorded  at  station  20S.  The  spectra  of  an  equal¬ 
ly  long  noise  window  prior  to  the  onset  of  P  are  also  shown.  The  am¬ 
plitude  values  are  in  arbitrary  units. 

11  Spectral  ratio  Salmon/Sterling,  corrected  for  noise,  with  logarithmic  18 

frequency  scale,  for  (a)  P  window  2.05  sec  long,  and  (b)  Lg  window 

4.10  sec  long,  recorded  at  station  20S.  The  dashed  lines  show  mean 
(least  squares)  slopes  over  the  frequency  range  5  to  25  Hz.  The  mean 
value  of  slope  and  its  standard  deviation  (S.D.)  are  indicted  for  fre¬ 
quency  ranges  of  5  to  25  and  5  to  20  Hz. 

12  Spectral  ratio  P/Lg,  corrected  for  noise,  for  (a)  Salmon  and  (b)  Sterling  19 

recorded  at  station  20S.  Points  for  which  S/N  power  ratio  is  less  than 

2  are  not  plotted.  The  dashed  lines  show  mean  (least  squares)  slopes 
over  the  frequency  range  5  to  25  Hz.  The  mean  value  of  slope/Hz  and 
its  standard  deviation  (S.D.)  are  indicted  for  frequency  ranges  of  5  to 
25  Hz  and  5  to  20  Hz. 

13  Spectral  ratio  P/Lg,  corrected  for  noise,  for  (a)  Salmon  and  (b)  Sterling  21 

recorded  at  station  PL-MS.  Points  for  which  S/N  power  ratio  is  less 

than  2  are  not  plotted.  The  dashed  lines  show  mean  (least  squares) 
slopes  over  the  frequency  range  5  to  25  Hz.  The  mean  value  of 
slope/Hz  and  its  standard  deviation  (S.D.)  are  indicated  for  frequency 
ranges  of  5  to  25  and  5  to  40  Hz. 

14  Spectra  of  the  observed  Pn  arrivals  (window  6.4  sec  long)  for  (a)  25 

Duryea  and  (b)  Buteo  recorded  at  station  KN-UT.  The  spectra  of  an 
equally  long  noise  window  prior  to  the  onset  of  Pn  are  also  shown. 

The  peak  amplitude  values,  noted  under  each  trace,  are  in  nanometers 
(nm). 
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15  Spectral  ratio  Duryea/Buteo,  corrected  for  noise,  with  logarithmic  fre-  26 
quency  scale,  for  (a)  Pn  window  6.4  sec  long,  and  (b)  Lg  window  12.8 

sec  long.  The  dashed  lines  show  mean  (least  squares)  slopes  over  the 
frequency  range  1  to  7  Hz.  The  mean  value  of  slope/Hz  and  its  stan¬ 
dard  deviation  (S.D.)  are  indicated  for  frequency  ranges  of  1  to  7  and  1 
to  5  Hz. 

16  Spectral  ratio  Pn/Lg,  corrected  for  noise,  for  (a)  Duryea  and  (b)  Buteo  27 

recorded  at  station  KN-UT.  Points  for  which  S/N  power  ratio  is  less 

than  2  are  not  plotted.  The  dashed  lines  show  mean  (least  squares) 
slopes  over  the  frequency  range  1  to  7  Hz.  The  mean  value  of 
slope/Hz  and  its  standard  deviation  (S.D.)  are  indicated  for  frequency 
ranges  of  1  to  7  Hz  and  1  to  5  Hz. 

17  Mean  Pn/Lg  spectral  slope  (per  Hz)  over  the  frequency  range  1  to  5  28 

Hz  versus  log  (scaled  depth)  for  7  explosions  (numbered  as  in  Table  1) 
recorded  at  station  KN-UT.  The  least  squares  linear  trend  (dashed 

line)  has  a  correlation  coefficient  of  0.810  and  slope  of  0.135.  The 
vertical  lines  represent  error  bars  with  one  standard  deviation. 

18  Spectra  corrected  for  noise,  source  function,  and  t*  for  (a)  Pn,  (b)  Lg,  30 

both  from  Duryea,  (c)  Pn,  and  (d)  Lg,  both  from  Buteo.  The  dashed 

lines  show  mean  (least  squares)  spectral  slopes  over  the  frequency 
range  1  to  5  Hz.  The  mean  value  of  slope/Hz  and  its  standard  devia¬ 
tion  (S.D.)  are  indicated  for  frequency  ranges  of  1  to  5  Hz  and  1  to  7 
Hz. 

19  Mean  spectral  slope  (per  Hz)  over  the  frequency  range  1  to  5  Hz  plot-  31 

ted  versus  log  (scaled  depth)  for  7  explosions  (numbered  as  in  Table  1) 
recorded  at  station  KN-UT  for  (a)  Pn,  the  least  squares  linear  trend 
(dashed  line)  has  correlation  coefficient  of  0.970  and  slope  of  0.206; 

and  (b)  Lg,  the  least  squares  linear  trend  (dashed  line)  has  correlation 
coefficient  of  0.662  and  slope  of  0.045. 

20  (a)  RDP’s  derived  from  the  observed  radial  horizontal  particle  velocity  35 

(URH)  at  E6-27  and  radial  particle  velocity  (UR)  at  E14C-36  with 
take-off  angles  of  90  and  31  degrees,  respectively.  Note  the  nearly 

equal  shotpoint-to-receiver  distances  indicated  for  the  two  receivers. 

(b)  RVP  spectral  ratio  E6-27/E14C-36,  corrected  for  noise  and  based 
on  only  about  0.2  sec  of  initial  ground  motion.  The  dashed  line  shows 
the  mean  (least  squares)  slope  of  0.00273  per  Hz  (with  standard  devia¬ 
tion  of  0.00014)  over  the  frequency  range  of  1  to  100  Hz.  The  inter¬ 
cept  value  is  1.032. 
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Most  specialists  on  seismic  verification  are  in  agreement  that  the  detonation  of  nuclear 
explosions  in  very  large  cavities  represents  the  most  significant  evasion  possibility  (Sykes, 
1987).  Cavity  decoupling  will  significantly  challenge  Comprehensive  Test  Ban  Treaty  (CTBT) 
verification  of  explosions  near  1  kt,  even  with  30  in-country  array  stations  (Hannon,  1986). 
The  largest  seismic  arrivals  recorded  by  these  arrays  will  be  the  crustal  wave  guide  phases 
such  as  Lg.  It  is  therefore  important  to  understand  the  differences  in  the  generation  of  P  and 
Lg  by  decoupled  and  tamped  explosions.  In  their  analysis  of  Salmon  and  Sterling  data,  Bland- 
ford  and  Woolson  (1979)  observed  decoupling  to  vary  as  a  function  of  frequency  and  the  fre¬ 
quency  dependence  was  substantially  different  for  P  and  Lg.  Their  findings  were,  however, 
essentially  based  on  data  from  a  single  station.  One  principal  objective  of  this  investigation  was 
to  verify  their  results  by  analyzing  data  from  other  stations  and  possibly  to  determine  why 
decoupling  varied  with  wave  type. 

The  decoupled  nuclear  explosion  Sterling  (December  3,  1966)  with  yield  of  0.38  kt  and 
shot  depth  of  829  m  was  shot  in  the  cavity  of  the  nuclear  explosion,  Salmon  (October  22, 
1964)  with  a  yield  of  5.3  kt  and  shot  depth  of  828  m.  Sterling  was  preceded  by  a  high- 
explosive  (HE)  calibration  shot  (November  17,  1966),  referred  to  hereafter  as  "Stercal",  with 
yield  equivalent  to  2.7  tons  of  TNT,  about  348  m  from  the  Sterling  location  and  at  about  the 
same  depth,  831  m  (Springer  et  al.,  1968).  The  three  explosions  were  recorded  at  several  sta¬ 
tions  at  varying  distances  and  azimuths  (see  Figure  1  and  Table  1  in  Springer  et  al.,  1968).  The 
instrument  responses  as  well  as  several  examples  of  recordings  can  be  seen  in  Springer  et  al. 
(1968)  and  Blandford  and  Woolson  (1979).  In  this  part  of  the  study,  we  analyzed  data  pri¬ 
marily  from  the  USCGS  stations  10S  and  20S  (16  and  32  km  south  of  shot  point)  since  data 
from  all  other  stations  had  already  been  examined  in  earlier  studies.  We  used  the  spectral  ratio 


Teledyne  Geotech 


I 


November  1986 


ANNUAL  REPORT 


STUDIES  IN  DECOUPLING 


TGAL-86-8 


method  in  which  the  spectra  of  observed  signals  from  various  sources,  after  correction  for 
noise,  are  compared  at  common  stations.  Care  was  taken  to  use  only  those  data  that  were  free 
from  spikes  as  well  as  clipped  signals. 
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ANALYSIS  OF  DATA  FROM  STATION  10S 

We  first  analyzed  the  Salmon,  Sterling,  and  Stercal  data  at  the  USCGS  station  10S 
located  about  10  miles  or  16  km  south  of  ground-zero.  Three-component  broadband  record¬ 
ings  of  the  three  explosions  at  one  of  several  receiver  positions  at  this  location  can  be  seen  in 
Springer  et  al.  (1968,  Figure  5).  Among  the  digitized  data  available  to  us,  at  common  sensors, 
only  the  vertical  component  data  were  found  to  be  usable.  The  best  quality  data  with  high  sig¬ 
nal  (S)  to  noise  (N)  ratio  at  station  10S  for  the  three  explosions  were  obtained  at  the  receiver 
location  C  2000;  these  are  shown  in  Figure  1  in  which  the  first  P  motions  are  aligned.  On 
each  record,  the  P  waves  are  followed  by  larger-amplitude,  lower  frequency  waves,  starting 
about  1.7  sec  after  the  onset  of  P.  On  the  basis  of  available  three-component  data  from  Sal¬ 
mon,  Sterling,  and  Stercal  at  comparable  source-receiver  distances  in  Borcherdt,  et  al.  (1967), 
this  prominent  arrival  is  too  early  to  be  the  direct  S  wave.  In  fact,  synthetic  seismograms  based 
on  source-receiver  distance,  A  =  16  km  and  crustal  structure  appropriate  for  this  region,  indi¬ 
cate  this  arrival  to  be  due  to  mostly  P  waves  reverberating  in  the  low-velocity  sediments  and 
the  S  waves  to  arrive  about  3.5  sec  after  P  (Charles  Langston,  written  communication).  The 
records  in  Figure  1  indicate  a  prominent  arrival  with  its  beginning,  denoted  by  S,  about  3.3  sec 
after  the  first  P  and  prominent  on  all  three  records.  Following  Blandford  and  Woolson  (1969), 
we  designate  this  group  of  waves  as  Lg.  In  spite  of  repeated  attempts,  the  data  could  not  be 
calibrated  correctly;  however,  as  we  shall  see  later,  this  did  not  adversely  affect  most  results  in 
this  study. 

The  term  Lg  is  generally  used  to  describe  the  prominent  regional  phase  consisting  of 
higher-mode  surface  waves  (Knopoff  et  al.,  1973).  Equivalently,  Lg  may  also  be  considered  to 
be  guided  shear  waves  incident  on  the  Moho  at  angles  more  grazing  than  the  critical  incidence 
and  multiply  reflected  within  the  crust  (Bouchon,  1982).  It  is  therefore  not  proper  to  refer  to 
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Figure  1.  Raw  common-receiver  vertical  component  data  recorded  at  station  10S  (Receiver 
location  C2000)  for  the  three  explosions  Salmon,  Sterling,  and  Stercal.  The  beginnings  of  the 
P  and  Lg  wavetrains,  indicated  by  the  letters  P  and  S,  are  separated  by  about  3.3  sec. 
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such  waves  as  Lg  at  source-receiver  distances  of  a  few  tens  of  km.  The  crustal  structure  in  the 
vicinity  of  Salmon  is  however  characterized  by  3.1  to  3.7  km  thick  coastal  sediments  (average 
P-wave  velocity  3  km/sec)  overlying  a  layer  with  velocity  of  5  km/sec  associated  with  strong 
refraction  arrivals  (Warren  et  al.,  1966).  This  provides  sufficient  impedance  contrast  for  the 
upper  low  velocity  layer  to  act  as  an  efficient  waveguide,  giving  rise  to  large-amplitude  waves 
in  the  same  manner  as  Lg.  The  numerical  modeling  of  Lg  by  Campillo  et  al.  (1985,  see  their 
Figures  9  and  10)  also  show  that,  for  a  shallow  source,  rays  which  contribute  to  the  Lg 
wavetrain  include  those  that  are  supercritically  reflected  on  the  upper,  shallower  discontinuities 
of  the  crust.  We  shall  therefore  retain  the  use  of  the  term  "Lg"  to  describe  the  low-frequency, 
large-amplitude  waves  such  as  those  observed  at  the  station  10S. 

The  digitization  rate  for  the  data  in  Figure  1  was  125  samples  per  sec.  For  each  record, 
the  initial  P  window  is  taken  to  be  about  1  sec  long  (actually  128  points),  a  10%  cosine  taper 
applied  and  the  amplitude  spectra  are  obtained  over  the  frequency  range  of  0  to  50  Hz.  This 
procedure,  the  same  as  used  by  Blandford  and  Woolson  (1979)  in  their  analysis  of 
Salmon/Sterling  data,  was  adopted  so  that  the  results  in  the  two  studies  can  be  directly  com¬ 
pared.  A  similar  procedure  was  used  for  obtaining  the  spectra  of  Lg  windows  which  were 
about  2  sec  long  (256  points);  the  starting  points  on  the  P  and  Lg  windows  were  separated  by 
about  1.6  sec.  These  spectra,  with  3  and  5  point  smoothing  for  P  and  Lg,  respectively,  are 
shown  in  Figures  2,  3,  and  4  for  Salmon,  Sterling,  and  Stercal,  respectively.  These  figures 
include  similar  spectra  of  samples  of  noise  prior  to  the  onset  of  P.  The  S/N  ratio  for  Salmon 
appears  to  be  good  only  up  to  about  25  Hz  whereas  for  Sterling  and  Stercal  S/N  is  good  for 
somewhat  higher  frequencies. 

The  decoupling  factor  is  defined  as  the  ratio  of  the  P-wave  amplitude  of  a  tamped  explo¬ 
sion  to  that  of  the  decoupled  shot  when  the  two  explosions  have  the  same  yield.  Since  Salmon 
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Spectra  of  the  observed  (a)  P  (window  1.02  sec  long)  and  (b)  Lg  (window  2.05  sec 
Salmon  recorded  at  10S.  The  spectra  of  noise  preceding  the  onset  of  direct  P  are 
n.  The  amplitude  values  are  in  arbitrary  units. 
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and  Sterling  had  quite  different  yields,  one  needs  to  divide  the  observed  Salmon/Sterling  ratio 
at  a  given  frequency  by  the  theoretical  spectral  ratio  of  a  5.3  kt  (Salmon)  explosion  to  a  0.38 
kt  (Sterling  tamped)  explosion.  We  obtained  the  theoretical  spectral  ratios  by  using  the  von 
Seggem  and  Blandford  (1972)  source  function  for  salt  (i.e.  k  =  26  at  5  kt  and  B  =  0).  The 
source  functions  for  explosions  with  yields  of  5.3,  0.38,  and  0.0027  kt,  shown  in  Figure  5, 
have  corner  frequencies  (half-amplitude  points)  at  about  4,  10,  and  50  Hz,  respectively.  The 
corresponding  spectral  ratios,  shown  in  Figure  6,  indicate  considerable  variation  with  fre¬ 
quency.  The  Salmon/Sterling  ratio  has  values  of  approximately  14  and  3  at  frequency  of  1 
and  20  Hz,  respectively,  whereas  Sterling/Stercal  ratio  is  approximately  140  and  30  at  frequen¬ 
cies  of  1  and  20  Hz,  respectively.  This  means  that,  over  the  frequency  range  of  1  to  20  Hz, 
each  of  the  theoretical  ratios  Salmon/Sterling(tamped)  and  Sterling(tamped)/Stercal  decreases 
by  a  factor  of  14/3  =  5.  Note  that  Springer  et  al.’s  (1968,  see  especially  their  Figure  17) 
theoretical  scaling  of  Salmon  and  Sterling  spectra  would  indicate  nearly  the  same  results  as  in 
Figure  6. 

The  spectral  ratios  Salmon/Sterling,  based  on  P  and  Lg  windows  containing  128  (1.02 
sec)  and  256  points  (2.05  sec),  respectively,  for  each  shot  and  recorded  at  the  same  common 
vertical  component  location  C  2000,  are  shown  in  Figures  7a  and  7b,  respectively.  The  ampli¬ 
tude  ratio  scales  are  obviously  in  error  (because  the  data  were  uncalibrated)  but  an  approxi¬ 
mately  correct  scale  was  derived  by  comparing  the  station  10S  amplitudes  in  Springer  et  al. 
(1968)  with  those  in  the  digitized  data;  these  are  indicated  at  the  left.  The  spectral  ratios  are 
corrected  for  noise  and  only  those  points  for  which  the  S/N  power  ratio  is  greater  than  2  for 
each  signal  are  included  and  the  dashed  line  shows  the  mean  (least  squares)  slope  over  the  fre¬ 
quency  range  of  5.0  to  25.0  Hz.  The  spectral  ratio  for  P  drops  at  a  fall-off  rate  of  about  f  _3  74 
over  the  frequency  range  of  5  to  25  Hz.  This  fall-off  rate  is  considerably  faster  than  about  f  ~2 
expected  on  the  basis  of  the  source  functions  in  Figure  6.  The  spectral  ratio  for  Lg,  shown  in 
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Figure  5.  Theoretical  spectra,  based  on  von  Seggern  and  Blandford’s  (1972)  source  function 
for  salt,  for  yields  of  5.3,  0.38,  and  0.0027  kt,  respectively.  The  corner  frequencies  (half 
amplitude  points)  for  the  three  spectra  are  at  about  4,  10,  and  50  Hz,  respectively. 
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Figure  6  Theoretical  spectral  ratios  for  Salmon/Sterling  (tamped),  Sterling  (tamped)/Stercal, 
and  Salmon/Stercal  over  the  frequency  range  1  to  100  Hz. 
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Figure  7.  Spectral  ratio  Salmon/Sterling,  corrected  for  noise,  with  logarithmic  frequency  scale, 
for  (a)  P  window  1.02  sec  long,  and  (b)  Lg  window  2.05  sec  long.  The  dashed  lines  show 
mean  (least  squares)  slopes  over  the  frequency  range  5  to  25  Hz.  The  mean  value  of  slope  and 
its  standard  deviation  (S.D.)  are  indicated. 
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Figure  7b,  indicate  a  fall-off  rate  of  about  f  ~2  94,  considerably  smaller  than  that  for  P,  over  the 
frequency  range  5  to  25  Hz.  The  P  amplitude  ratio  (Figure  7a)  drops  by  a  factor  of  about 
2000/30  =  70  between  5  and  20  Hz.  Similarly,  in  Figure  7b,  the  Lg  amplitude  ratio  drops  by  a 
factor  of  about  2000/50  =  40  between  5  and  20  Hz.  Thus,  at  higher  frequencies,  the  Sterling 
Lg  is  more  decoupled  than  the  Sterling  P  wave.  The  amplitude  ratios  at  low  frequencies  and 
the  significantly  different  fall-off  rates  for  P  and  Lg  are  similar  to  those  obtained  by  Blandford 
and  Woolson  (1979,  see  their  Figure  15). 

The  observed  spectral  ratios  for  P  in  Figure  7(a)  suggest  the  Salmon/Sterling  ratios  to  be 
about  2000  and  30  at  frequencies  of  1  Hz  (i.e.  low-frequency  value)  and  20  Hz,  respectively. 
The  corrected  decoupling  factors  for  Salmon/Sterling  are  therefore  2000/14  =  140  and  30/3  = 
10,  at  frequencies  of  1  and  20  Hz,  respectively.  The  Sterling  decoupling  ratio  for  P  is  there¬ 
fore  reduced  by  a  factor  of  14  over  the  frequency  range  of  1  to  20  Hz.  Similarly,  on  the  basis 
of  Figure  7b,  the  corrected  decoupling  factors  for  Lg  are  approximately  140  and  17  at  frequen¬ 
cies  of  1  and  20  Hz,  respectively,  and  the  decoupling  ratio  is  reduced  by  a  factor  of  only  about 
8  over  the  frequency  range  of  1  to  20  Hz. 

The  spectral  ratios  P/Lg  for  the  three  explosions,  again  from  the  same  common  vertical 
component  location  and  corrected  for  noise,  are  shown  in  Figure  8.  Note  the  similar  spectral 
modulations  in  the  three  plots  for  frequencies  up  to  about  20  Hz;  this  is  probably  because  the 
source-receiver  paths  are  almost  identical  in  the  three  cases.  Comparing  the  results  for  Salmon 
and  Sterling,  the  differences  in  the  two  spectral  ratios  exist  primarily  at  the  higher  frequencies. 
The  dashed  lines  in  Figure  8  show  the  mean  (least  squares)  slopes  over  the  frequency  range  of 
5  to  25  Hz.  For  each  explosion,  the  mean  spectral  slope  values  with  its  standard  deviation  (S. 
D.)  are  indicated  for  frequency  bands  of  5  to  25  Hz  and  5  to  40  Hz.  The  mean  slope  values 
are  significantly  different  for  Salmon  and  Sterling  and  are  the  smallest  for  Salmon  and  the 
largest  for  Stercal. 
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Figure  8.  Spectral  ratio  P/Lg,  corrected  for  noise,  for  (a)  Salmon,  (b)  Sterling  and  (c)  Stercal 
recorded  at  the  station  10S.  Points  for  which  S/N  power  ratio  is  less  than  2  are  not  plotted. 
The  dashed  lines  show  mean  (least  squares)  slopes  over  the  frequency  range  5  to  25  Hz.  The 
mean  value  of  slope/Hz  and  its  standard  deviation  (S.D.)  are  indicated  for  frequency  ranges  of 
5  to  25  Hz  and  5  to  40  Hz. 
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ANALYSIS  OF  DATA  FROM  STATION  20S 

Recordings  at  the  station  20S,  located  about  20  miles  or  32  km  south-southwest  of  the 
ground-zero,  provided  vertical  component  data  with  good  S/N  for  both  Salmon  and  Sterling. 
Nearly  24  sec  of  common-receiver  data  are  shown  in  Figure  9  in  which,  similar  to  Figure  1, 
the  P  waves  are  followed  by  larger-amplitude,  lower-frequency  wave-packet,  Lg,  the  onset  of 
which  is  denoted  by  S.  The  beginnings  of  P  and  Lg  are  separated  by  about  5.5  sec  which 
agrees  well  with  the  synthetics  for  A  =  32  km  (Charles  Langston,  written  communication).  The 
spectra  of  P  and  noise  windows,  each  256  points  (2.05  sec)  long,  for  Salmon  and  Sterling, 
obtained  in  the  same  manner  as  for  the  data  from  station  10S,  are  shown  in  Figure  10.  The 
spectral  ratios  Salmon/Sterling  for  P  and  Lg  windows  of  lengths  256  and  512  points,  respec¬ 
tively,  are  shown  in  Figure  11.  Over  the  frequency  range  of  5  to  25  Hz,  the  amplitude  ratios 
have  fall-off  rates  of  about  f-301  and  f-1'44  for  P  and  Lg,  respectively.  The  higher  fall-off  rate 
for  the  P  ratio  than  for  the  Lg  ratio  is  similar  to  the  results  obtained  from  the  data  at  station 
10S  (Figure  7). 

The  spectral  ratios  P/Lg  for  Salmon  and  Sterling,  corrected  for  noise,  are  shown  in  Figure 
12.  The  spectral  modulations  for  the  two  explosions  are  similar,  especially  at  the  lower  fre¬ 
quencies.  The  differences  in  the  spectral  ratios  lie  mainly  in  their  higher  frequency  values.  The 
dashed  lines,  representing  mean  slope  values  over  the  frequency  range  of  5  to  25  Hz,  show 
considerably  smaller  (more  negative)  slope  for  Salmon  than  for  Sterling;  the  difference  between 
the  two  slope  values  is  0.0414.  All  these  results  are  similar  to  those  for  data  from  station  10S 
(Figure  8). 
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Figure  9.  Raw  common-receiver,  vertical  component  data  at  station  20S  for  Salmon  and  Ster¬ 
ling.  The  beginnings  of  the  P  and  Lg  wavetrains,  denoted  by  the  letters  P  and  S,  are  separated 
by  about  5.5  sec. 
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Figure  10.  Spectra  of  the  observed  P  arrivals  (window  2.05  sec  long)  for  (a)  Salmon  and  (d) 
Sterling  recorded  at  station  20S.  The  spectra  of  an  equally  long  noise  window  prior  to  the 
onset  of  P  are  also  shown.  The  amplitude  values  are  in  arbitrary  units. 
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figure  11.  Spectral  ratio  Salmon/Sterling,  corrected  for  noise,  with  logarithmic  frequency 
scale,  for  (a)  P  window  2.05  sec  long,  and  (b)  Lg  window  4.1U  sec  long,  recorded  at  station 
20S.  The  dashed  lines  show  mean  (least  squares)  slopes  over  the  frequency  range  5  to  25  Hz. 
The  mean  value  of  slope  and  its  standard  deviation  (S.D.)  are  indicted  for  frequency  ranges  of 
5  to  25  and  5  to  20  Hz. 
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Figure  12.  Spectral  ratio  P/Lg,  corrected  for  noise,  for  (a)  Salmon  and  (b)  Sterling  recorded  at 
station  20S.  Points  for  which  S/N  power  ratio  is  less  than  2  are  not  plotted.  The  dashed  lines 
show  mean  (least  squares)  slopes  over  the  frequency  range  5  to  25  Hz.  The  mean  value  of 
slope/Hz  and  its  standard  deviation  (S.D.)  are  indicted  for  frequency  ranges  of  5  to  25  Hz  and 
5  to  20  Hz. 
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ANALYSIS  OF  DATA  FROM  STATION  PL-MS 

We  also  analyzed  the  vertical  component  data  for  Salmon  and  Sterling  recorded  at  the 
USGS  station  PL-MS,  located  27  km  south  of  ground-zero.  The  digitization  rate  for  the  PL- 
MS  data  was  100  samples/sec.  The  spectral  ratios  P/Lg  for  Salmon  and  Sterling,  with  the  Lg 
window  being  512  points  (5.12  sec)  long,  are  shown  in  Figure  13.  The  spectral  modulations 
are  again  similar,  especially  for  the  lower  frequencies.  The  values  of  mean  slopes  (dashed 
lines)  are  significantly  smaller  for  Salmon  than  for  Sterling.  The  difference  between  the  two 
mean  slope  values  is  0.0463  for  the  frequency  range  of  5  to  25  Hz,  nearly  equal  to  the 
difference  observed  at  station  20S.  These  results  are  again  similar  to  those  obtained  from  data 
at  stations  10S  and  20S. 
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Figure  13.  Spectral  ratio  P/Lg,  corrected  for  noise,  for  (a)  Salmon  and  (b)  Sterling  recorded  at 
station  PL-MS.  Points  for  which  S/N  power  ratio  is  less  than  2  are  not  plotted.  The  dashed 
lines  show  mean  (least  squares)  slopes  over  the  frequency  range  5  to  25  Hz.  The  mean  value 
of  slope/Hz  and  its  standard  deviation  (S.D.)  are  indicated  for  frequency  ranges  of  5  to  25  and 
5  to  40  Hz. 


y 


ANNUAL  REPORT  STUDIES  IN  DECOUPLING  TGAL-86-8 

SPECTRAL  ANALYSIS  OF  Pn  AND  Lg  FROM  NTS  EXPLOSIONS 

An  examination  of  P  and  the  Lg  group  data  from  the  three  stations  10S,  20S,  and  PL- 
MS,  described  above,  shows  exactly  similar  differences  in  the  excitation  of  the  two  phases 
from  Salmon  and  Sterling.  An  obvious  major  difference  between  Salmon  and  Sterling  is  in 
their  scaled  depths,  defined  as  depth(m)AVu3  where  W  is  yield  in  kt.  In  order  to  understand 
why  the  decoupling  ratio  is  different  for  P  and  Lg,  we  therefore  studied  the  spectral  content  of 
Pn  and  Lg  from  explosions  covering  a  wide  range  of  yields,  shot  depths,  and  scaled  depths. 
For  this  purpose,  we  selected  7  Pahute  Mesa  (NTS)  explosions,  covering  the  yield  range  of 
0.7  to  1100  kt  (Table  1).  Regional  phases  from  these  explosions  were  well  recorded  at  the 
LRSM  station,  KN-UT  (epicentral  distance,  A  =  320  km);  Pn  arrivals  from  these  were  studied 
by  Blandford  (1976).  Two  of  the  seven  explosions,  viz.  Cabriolet  and  Schooner,  were  cratering 
shots.  It  is  interesting  to  note  that  time-domain  measurements  of  the  amplitudes  of  regional 
phases  (Pn,  Pg,  and  Lg)  from  these  7  explosions  failed  to  show  any  systematic  differences 
between  the  shallowest  (cratering)  and  the  deepest  and/or  the  most  overburied  explosions 
(Gupta  et  al.,  1985). 
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TABLE  1 


PAHUTE  MESA  AND  SALT  EXPLOSIONS  USED  IN  SPECTRAL  ANALYSES 


No.  Name  Date  Shot  Shot  Depth  of  Yield+  Scaled 

Medium  Depth  Water  Table  Depth# 

(m)  (m)  (kt) 

1 

Buteo 

12  May  65 

i 

tuff 

696 

660 

0.7 

784 

2 

Rex 

24  Feb  66 

tuff 

671 

642 

16 

266 

3 

Duryea 

14  Apr  66 

rhyolite 

544 

662 

65 

135 

4 

Scotch 

23  May  67 

tuff 

977 

672 

150 

184 

5 

Cabriolet* 

26  Jan  68 

rhyolite 

52 

488 

2.3 

39 

6 

Schooner* 

8  Dec  68 

tuff 

111 

274 

35 

34 

7 

Benham 

19  Dec  68 

tuff 

1402 

641 

1100 

136 

8 

Salmon 

22  Oct  64 

salt 

828 

- 

5.3 

475 

9 

Stercal 

17  Nov  66 

salt 

831 

- 

0.0027 

5968 

10 

Sterling 

3  Dec  66 

salt 

829 

- 

0.38 

1145 

+  Values  for  NTS  explosions  from  Springer  and  Kinnaman  (1971)  except  for  Buteo  from  an 
estimate,  based  on  LR  data,  by  Blandford  (1976).  Values  for  the  three  shots  in  salt  from 
Springer  et  al.  (1968). 

*  cratering  explosion 

#  depth(m)/W1/3  where  W  is  yield  in  kt. 

Scaled  depth  is  an  important  parameter  directly  influencing  the  spectra  and  amplitudes  of 
explosion-generated  body  and  surface  waves  (e.g.  Rulev,  1965;  and  Kisslinger,  1963,  for  a  re¬ 
view  of  early  work).  Table  1  lists  the  scaled  depths  not  only  for  the  7  Pahute  Mesa  explosions 
but  also  for  Salmon,  Sterling,  and  Stercal.  A  large  range  of  both  shot  and  scaled  depths  are 
covered  by  these  10  explosions. 

We  obtained  the  spectral  ratios  of  Pn/Lg  on  the  vertical  component  records  at  the  station, 
KN-UT  for  the  seven  Pahute  Mesa  (NTS)  explosions  listed  in  Table  1.  To  isolate  the  source 
effects,  it  is  important  to  use  as  nearly  common  source-receiver  paths  as  possible  because  of 
the  extreme  sensitivity  to  propagation  paths  of  the  regional  phases,  especially  Lg  (Gupta  and 
Blandford,  1983).  The  seven  explosions  in  Table  1  are  all  from  Pahute  Mesa  so  that  the 
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source-receiver  paths  to  KN-UT  should  not  be  much  different  and  the  results  should  be  nearly 
free  from  propagation  path  effects.  The  window  lengths  for  Pn  and  Lg  were  6.4  and  12.8  sec, 
respectively,  and  a  10%  cosine  taper  was  used  along  with  3  and  5  point  smoothing  of  the  spec¬ 
tra.  Out  of  these  7  shots,  Buteo  and  Duryea  with  their  yields  differing  by  a  factor  of  about  100, 
were  detonated  in  the  same  hole  so  that  the  source-receiver  paths  are  almost  identical.  The 
spectra  of  Pn  and  noise  windows  for  Duryea  and  Buteo  are  shown  in  Figure  14  which  show 
good  S/N  up  to  about  7  Hz.  The  spectral  ratios  Duryea/Buteo  for  Pn  and  Lg  windows,  each 
corrected  for  noise,  are  shown  in  Figure  15.  The  fall-off  with  frequency  rate  is  larger  for  the 
Pn  ratio  than  for  the  Lg  ratio,  similar  to  the  results  for  Salmon/Sterling  in  Figures  7  and  11. 
The  spectral  ratios  Pn/Lg  for  Duryea  and  Buteo,  corrected  for  noise,  are  shown  in  Figure  16 
which  indicates  smaller  spectral  slope  for  Duryea  than  for  Buteo,  again  similar  to  the  Salmon 
and  Sterling  results  of  Figures  8,  12,  and  13.  Results  from  all  7  explosions  (Table  1)  are 
shown  in  Figure  17  in  which  the  mean  spectral  slopes  of  Pn/Lg  over  the  frequency  range  of  1 
to  5  Hz  are  plotted  versus  the  scaled  depth  (in  log  units).  The  plot  indicates  a  linear  trend  with 
correlation  coefficient  of  0.810.  This  implies  that  the  spectral  ratio  Pn/Lg  is  considerably 
influenced  by  scaled  depth. 

It  is  important  to  determine  whether  it  is  the  spectra  of  Pn  or  Lg  (or  both)  that  vary 
strongly  with  scaled  depth.  A  comparison  of  the  spectral  ratios  for  P  or  Pn  and  Lg,  such  as  in 
Figures  7  and  11  for  Salmon/Sterling  and  in  Figure  15  for  Duryea/Buteo,  suggests  that  the  Lg 
spectra  are  less  variable  or  more  stable  than  the  P  or  Pn  spectra.  To  explore  this  possibility  on 
the  observed  data  from  the  7  Pahute  Mesa  explosions  recorded  at  KN-UT,  we  compared  the 
spectral  content  of  the  individual  Pn  and  Lg  windows  after  correcting  for  the  large  differences 
in  the  explosion  yields.  For  the  latter,  we  made  use  of  von  Seggern  and  Blandford’s  (1972) 
source  scaling  in  which  the  source  function  for  a  given  yield,  W  is  expressed  as 
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Figure  14.  Spectra  of  the  observed  Pn  arrivals  (window  6.4  sec  long)  for  (a)  Duryea  and  (b) 
Buteo  recorded  at  station  KN-UT.  The  spectra  of  an  equally  long  noise  window  prior  to  the 
onset  of  Pn  are  also  shown.  The  peak  amplitude  values,  noted  under  each  trace,  are  in  nanom¬ 
eters  (nm). 
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Figure  15.  Spectral  ratio  Duryea/Buteo,  corrected  for  noise,  with  logarithmic  frequency  scale, 
for  (a)  Pn  window  6.4  sec  long,  and  (b)  Lg  window  12.8  sec  long.  The  dashed  lines  show 
mean  (least  squares)  slopes  over  the  frequency  range  1  to  7  Hz.  The  mean  value  of  slope/Hz 
and  its  standard  deviation  (S.D.)  are  indicated  for  frequency  ranges  of  1  to  7  and  1  to  5  Hz. 
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Figure  16.  Spectral  ratio  Pn/Lg,  corrected  tor  noise,  for  (a)  Duryea  and  (b)  Buteo  recorded  at 
station  KN-UT.  Points  for  which  S/N  power  ratio  is  less  than  2  are  not  plotted.  The  dashed 
lines  show  mean  (least  squares)  slopes  over  the  frequency  range  1  to  7  Hz.  The  mean  value  of 
slope/Hz  and  its  standard  deviation  (S.D.)  are  indicated  for  frequency  ranges  of  1  to  7  Hz  and 
1  to  5  Hz. 
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Figure  17.  Mean  Pn/Lg  spectral  slope  (per  Hz)  over  the  frequency  range  1  to  5  Hz  versus  log 
(scaled  depth)  for  7  explosions  (numbered  as  in  Table  1)  recorded  at  station  KN-UT.  The  least 
squares  linear  trend  (dashed  line)  has  a  correlation  coefficient  of  0.810  and  slope  of  0.135. 
The  vertical  lines  represent  error  bars  with  one  standard  deviation. 
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S(  co)  =  W 


A2((0/k)2  +  1 


(oVk)2  +  1 


1/2 

iir  ^(~> 


where  to  =  2  n  f,  'F(oo)  is  a  coupling  term  independent  of  frequency; 


k  =  ko  (5/W) 


1/3 


and 


(1) 


(2) 


A  =  2B  +  1  (3) 

where  B  and  kQ  are  medium-dependent  constants.  In  this  model,  shot  depth  is  not  an  indepen¬ 
dent  parameter  but  the  source  function  can  be  varied  by  assigning  different  values  to  the  rise 
time  and  overshoot  parameters  ko  and  B,  respectively.  Considering  first  the  Pn  windows  for  the 
various  explosions,  the  spectra  were  corrected  for  the  instrument  response  and  anelastic  at¬ 
tenuation  by  assuming  t*  =  0.1,  the  value  used  by  Blanaford  (1976).  A  correction  for  the 
differences  in  yields  was  applied  by  using  the  source  function  in  equation  (1)  corresponding  to 
a  tuff  reduced  displacement  potential  (RDP),  i.e.  ko  =  12.0  and  B  =  0.  With  these  three 
corrections,  the  Pn  spectra  would  become  flat  if  the  shot  depth  and  medium  did  not  influence 
the  spectral  shapes.  Most  of  the  corrected  spectra  were  not  flat  and  their  mean  spectral  slopes 
indicated  nearly  linear  trends;  the  results  for  Duryea  and  Buteo  are  shown  in  Figure  18.  This 
means  that  ordinary  cube-root  scaling  (i.e.  with  fixed  values  of  parameters  ko  and  B)  fails  to 
model  the  observed  Pn  spectra.  A  plot  of  the  derived  mean  spectral  slopes  (such  as  those  in 
Figure  18),  over  the  frequency  range  1.0  to  5.0  Hz,  versus  the  scaled  depth  (in  log  units)  is 
shown  in  Figure  19a.  The  least  squares  linear  trend  (dashed  line)  has  correlation  coefficient  of 
0.970  with  mean  slope  of  0.206  indicating  that,  for  a  fixed  yield,  the  Pn  spectrum  of  an  over¬ 
buried  (large  scaled  depth)  explosion  is  richer  in  higher  frequencies  compared  to  that  of  an  un¬ 
derburied  (small  scaled  depth)  explosion.  This  is  in  qualitative  agreement  with  Mueller  and 
Murphy’s  (1971)  scaling  in  which  scaled  depth  is  not  an  important  parameter  but  shot  depth  is 
relevant  because  of  its  influence  on  the  overburden  pressure.  The  same  corrections  (i.e.  t*  = 
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Figure  18.  Spectra  corrected  for  noise,  source  function,  and  t*  for  (a)  Pn,  (b)  Lg,  both  from 
Duryea,  (c)  Pn,  and  (d)  Lg,  both  from  Buteo.  The  dashed  lines  show  mean  (least  squares) 
spectral  slopes  over  the  frequency  range  1  to  5  Hz.  The  mean  value  of  slope/Hz  and  its  stan¬ 
dard  deviation  (S.D.)  are  indicated  for  frequency  ranges  of  1  to  5  Hz  and  1  to  7  Hz. 
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Figure  19.  Mean  spectral  slope  (per  Hz)  over  the  frequency  range  1  to  5  Hz  plotted  versus 
log  (scaled  depth)  for  7  explosions  (numbered  as  in  Table  1)  recorded  at  station  KN-UT  for  (a) 
Pn,  the  least  squares  linear  trend  (dashed  line)  has  correlation  coefficient  of  0.970  and  slope  of 
0.206;  and  (b)  Lg,  the  least  squares  linear  trend  (dashed  line)  has  correlation  coefficient  of 
0.662  and  slope  of  0.045. 
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0.1  and  tuff  RDP)  were  also  applied  to  the  Lg  spectra  and  the  mean  spectral  slopes  over  the 
same  frequency  range  computed.  A  plot  of  these  slopes  versus  the  scaled  depth  is  shown  in 
Figure  19b;  the  linear  trend  (dashed  line)  has  correlation  coefficient  of  0.662  and  mean  slope  of 
0.045.  Note  that  whereas  the  Pn  slope  values  in  Figure  19a  cover  a  range  of  about  0.28  per 
Hz,  the  Lg  slope  values  cover  a  range  of  only  about  0.08  per  Hz.  This,  together  with  the  mean 
slope  value  for  the  data  in  Figure  19a  being  considerably  larger  than  that  in  Figure  19b,  sug¬ 
gest  that,  for  a  fixed  yield,  the  Lg  spectra  (unlike  the  Pn  spectra),  are  insensitive  to  variations 
in  scaled  depth. 

Comparing  results  for  the  mean  spectral  slopes  of  Pn  and  Lg  separately  (Figures  19a  and 
19b)  with  those  for  Pn/Lg  (Figure  17),  it  appears  that  most  of  the  observed  spectral  differences 
between  Pn  and  Lg  are  due  to  the  spectra  of  P  waves  varying  significantly  with  scaled  depth 
whereas  the  spectra  of  Lg  are  insensitive  to  scaled  depth.  In  other  words,  the  observed  varia¬ 
tion  in  the  spectral  slope  Pn/Lg  is  essentially  due  to  the  P-wave  spectra  being  significantly 
more  dependent  on  the  scaled  depth  than  the  Lg  spectra. 
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POSSIBLE  EXPLANATION  FOR  VARIATION  IN  P/Lg  FOR  SALMON  AND  STERLING 

At  each  of  the  three  stations  10S,  20S,  and  PL-MS,  the  spectral  ratio  P/Lg  was  found  to 
increase  with  increasing  frequency,  more  positively  for  Sterling  than  for  Salmon,  i.e.  the  spec¬ 
tral  slope  was  larger  for  larger  scaled  depth,  in  agreement  with  the  results  from  7  Pahute  Mesa 
explosions  recorded  at  KN-UT.  Moreover,  a  comparison  of  the  spectral  ratios  Salmon/Sterling 
for  P  and  Lg  suggested  that  the  P  spectra  were  considerably  more  variable  (less  stable)  than 
the  Lg  spectra,  again  in  agreement  with  the  results  from  Pahute  Mesa  data.  A  possible  expla¬ 
nation  for  the  spectral  ratio  P/Lg  to  be  different  for  P  and  Lg  is  therefore  that  the  dependence 
on  scaled  depth  of  the  P  spectra  is  considerably  larger  than  for  the  Lg  group. 

The  scaled  depths  for  Sterling  and  Salmon  differ  by  a  factor  of  2.41  (see  Table  1).  If  we 
assume  that  P/Lg  =  Pn/Lg,  then,  on  the  basis  of  the  empirical  relationship  in  Figure  17,  the 
mean  spectral  slopes  P/Lg  for  the  two  explosions  should  differ  by  about  0.05/Hz.  The  ob¬ 
served  differences  from  data  at  stations  10S,  20S,  and  PL-MS,  over  the  frequency  range  of  5 
to  25  Hz,  are  0.015/Hz,  0.041/Hz,  and  0.046/Hz,  respectively  (see  Figures  8,  12,  and  13).  The 
agreement  between  expected  and  observed  values  is  good  for  stations  20S  and  PL-MS  but  poor 
for  station  10S,  probably  because  the  Lg  group  has  not  fully  developed  due  to  the  small 
source-receiver  distance.  Note  that  Lg  is  the  largest-amplitude  phase  for  data  at  stations  20S 
and  PL-MS  but  not  at  10S. 

A  possible  but  not  necessarily  valid  mechanism  for  the  variability  of  Pn  and  stability  of 
Lg  spectra  may  be  offered  on  the  basis  of  RDP  varying  significantly  with  the  takeoff  angle, 
proposed  earlier  by  Blandford  (1976).  The  explosion  cavity  associated  with  an  overburied  shot 
is  likely  to  be  nearly  spherical  whereas  that  for  a  shot  with  smaller  scaled  depth  is  expected  to 
approximate  an  ellipsoid  with  the  long  axis  along  the  vertical  direction.  In  terms  of  moment- 
tensor  source  (Aki  and  Richards,  1980,  Chapter  3),  a  spherical  cavity  has  M^M^Mjj 
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whereas  the  equalities  do  not  hold  for  non-spherical  cavity.  A  source  with  azimuthal  symmetry 
may  have  Mm  =  M22  *■  M33  so  that  a  cylindrically  symmetric  source  extended  in  the  vertical 
direction  would  appear,  in  the  far  field,  to  have  the  time  function  for  M33  different  than  for 
Mm  and  M22.  For  example,  the  corner  frequencies  and  the  DC  levels  of  Mn  and  M22  would 
be  higher  than  for  M33  for  a  pressurized  ellipsoid  with  major  axis  in  the  vertical  direction 
(Glenn  et  al.,  1986,  especially  their  Figure  5). 

Near- field  subsurface  data  from  nuclear  explosions,  sufficient  to  test  the  above  possibility, 
perhaps  do  not  exist.  Limited  data  from  Salmon,  based  on  just  two  detectors,  do  appear  to  sup¬ 
port  the  proposed  mechanism,  in  spite  of  Salmon  being  an  overburied  shot  with  small  expected 
source  anisotropy.  Figure  20a  (after  Figure  5.20,  Perret,  1968)  shows  the  measured  RDP’s 
based  on  the  observed  particle  velocities  at  two  detectors  E6-27  and  E14C-36,  with  takeoff  an¬ 
gles,  measured  from  the  vertical,  of  90  and  31  deg,  respectively.  The  source-receiver  distances 
are  nearly  equal  so  that  the  substantial  difference  in  the  two  RDP’s  cannot  be  attributed  to  spa¬ 
tial  attenuation.  The  initial  pulse  rises  somewhat  faster  for  the  shot-level  receiver,  E6-27  than 
for  the  other  receiver,  E14C-36.  Moreover,  the  static  RDP  value,  which  corresponds  to  the  DC 
level  in  the  far-field,  is  higher  for  E6-27  than  for  E14C-36.  Similar  results  are  obtained  by 
analyzing  the  reduced  velocity  potential  (RVP),  which  corresponds  to  the  observed  ground 
motion  in  the  far-field.  The  RVP  spectral  ratio  E6-27/E14C-36,  based  on  signal  windows  of 
only  about  0.2  sec  each,  is  shown  in  Figure  20b.  The  effect  of  the  free  surface  reflected  phases 
(see  Figure  10,  Glenn  et  al.,  1986)  should  be  minimal  since  the  signal  windows  are  much 
smaller  than  the  teleseismic  pP  -  P  time  of  about  0.58  sec.  In  Figure  20b,  the  source  time 
function  along  the  horizontal  direction  (E6-27)  is  richer  in  higher  frequencies  and  has  some¬ 
what  (about  3%)  larger  DC  amplitude  than  that  along  the  near- vertical  direction  (E14C-36),  in 
qualitative  agreement  with  Glenn  et  al.  (1986). 
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Pn  consists  mainly  of  waves  leaving  the  focal  sphere  at  steep  (nearly  vertical)  angles 
whereas  Lg  spans  a  large  range  of  takeoff  angles  around  the  horizontal  direction  (see  e.g.  Cam- 
pillo  et  al.,  1985).  Thus  Pn  is  sensitive  to  variations  in  the  M33  component  whereas  Lg  is  less 
sensitive  to  M33  because  it  averages  over  Mu,  M22,  and  M33.  For  shots  with  small  scaled 
depth,  the  effect  of  the  source  cavity  will  therefore  be  to  shift  the  spectra  of  Pn  to  lower  fre¬ 
quencies  and  that  of  Lg  to  relatively  higher  frequencies.  Another  contributing  factor  can  be  the 
non-linear  free-surface  effects,  more  intense  for  shots  with  smaller  scaled  depth  and  mostly 
confined  to  the  region  between  the  shot  point  and  the  free  surface  (e.  g.  Rimer  et  al.,  1979, 
especially  Figures  2.67  and  2.68).  Vertically  traveling  waves  that  contribute  directly  to  Pn, 
such  as  pP,  will  pass  through  this  highly-fractured,  large-strain,  low-Q  region  and  may  become 
deficient  in  high  frequencies.  Frasier  (1972)  and  Shumway  and  Blandford  (1978)  provided  evi¬ 
dence  that  pP  is  deficient  in  higher  frequencies  as  compared  to  direct  P.  Lg  should  be  relatively 
free  from  these  free-surface  effects.  We  can  thus  qualitatively  explain  the  observed  increase  in 
the  Pn/Lg  spectral  slope  with  scaled  depth  such  as  that  shown  in  Figure  17. 
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DISCUSSION  AND  CONCLUSIONS 

A  comparison  of  the  P/Lg  spectral  ratios  for  Salmon  and  Sterling  at  stations  10S  (Figure 
8),  20S  (Figure  12),  and  PL-MS  (Figure  13)  shows  that  most  of  the  differences  in  the  spectra 
of  P  and  Lg  lie  in  their  high  frequency  contents.  Similarly,  for  the  NTS  explosions  recorded 
at  KN-UT,  the  mean  spectral  slope  of  the  spectral  ratio  Pn/Lg  increases  with  scaled  depth 
mostly  because  of  differences  in  the  Pn  and  Lg  spectra  at  high  frequencies  (see  Figure  16).  In 
other  words,  the  low  frequencies  in  P  and  Lg  are  insensitive  to  variations  in  scaled  depth 
whereas  the  high  frequencies  vary  considerably  with  scaled  depth.  Flynn  (1986)  studied  effects 
of  source  depth  on  near  source  seismograms  by  examining  the  spectra  of  P  and  SV-Rayleigh 
surface  waves  from  five  115  kg  chemical  explosions  buried  at  depths  ranging  from  1.8m  (op¬ 
timum  cratering  depth)  to  11.5  m  (fully  contained).  Energy  in  the  high  frequency  band  of  10  to 
40  Hz,  dominated  by  P  waves,  was  found  to  increase  significantly  with  source  depth  whereas 
energy  in  the  low  frequency  band  of  4  to  12  Hz,  consisting  mostly  of  SV-Rayleigh  waves,  was 
not  a  strong  function  of  shot  depth.  Our  results  regarding  the  relative  excitation  of  P  and  Lg 
are  therefore  in  good  agreement  with  those  of  Flynn  (1986). 

Underground  nuclear  explosions  are  generally  very  complex  with  both  primary  and 
secondary  sources  of  seismic  waves  (e.g.  Masse,  1981;  Springer,  1974).  Langston  (1983) 
found  a  complex  mechanism  even  for  the  decoupled  shot  Sterling.  Blandford  and  Woolson 
(1979)  noted  evidence  for  possible  azimuthal  variation  in  the  decoupling  ratio  for  Salmon  and 
Sterling.  The  idealization  of  an  explosion  source  as  a  pressure  pulse  applied  to  the  walls  of  a 
spherical  cavity  should  perhaps  be  used  to  infer  only  qualitative  results.  Non-linear  effects, 
known  to  be  important  for  most  nuclear  explosions  (especially  those  with  shallower  scaled 
depths),  also  drastically  modify  the  seismic  signals  (Perl  et  al.,  1979;  Rimer  et  al.,  1979;  Day 
et  al.,  1985). 
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Analysis  of  data  from  Salmon,  Sterling,  and  the  Sterling  Calibration  shot  recorded  at  both 
10S  and  20S  indicate  reduced  decoupling  at  higher  frequencies  for  both  P  and  the  Lg  group. 
After  correcting  for  the  differences  in  yield  of  the  three  explosions,  the  reduction  in  decoupling 
for  P  is  a  factor  of  about  7  to  14  over  the  frequency  range  1  to  20  Hz.  The  reduction  in 
decoupling  for  Lg  over  the  same  frequency  range  is  considerably  smaller.  Therefore,  at  high 
frequencies,  Lg  is  more  decoupled  than  P.  These  results  are  in  agreement  with  the  earlier  work 
of  Blandford  and  Woolson(1979).  Note  that  the  absolute  calibration  of  data  used  in  this  study 
is  somewhat  uncertain  so  that  the  earlier  estimates  of  the  Salmon/Sterling  decoupling  ratio  by 
Blandford  and  Woolson  (1979)  are  to  be  preferred. 

A  comparison  of  the  spectra  of  the  direct  P  and  the  Lg  group  for  Salmon  and  Sterling 
recorded  at  each  of  the  three  stations  10S,  20S,  and  PL-MS  suggests  the  P-wave  spectra  to  be 
significantly  more  variable  than  the  Lg  spectra.  Analysis  of  Pn  and  Lg  spectra  from  several 
Pahute  Mesa  shots  recorded  at  a  common  station  indicates  similar  results  and  a  rather  sys¬ 
tematic  dependence  of  the  mean  slope  of  the  spectral  ratio  Pn/Lg  on  scaled  depth.  For  a  fixed 
yield,  Pn  spectra  vary  considerably  with  scaled  depth  whereas  the  Lg  spectra  are  insensitive  to 
scaled  depth.  Differences  in  the  frequency  dependence  of  decoupling  for  P  and  Lg  are  there¬ 
fore  probably  due  to  the  P  spectra  varying  strongly  with  scaled  depth  and  the  Lg  spectra  being 
nearly  independent  of  the  scaled  depth. 
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SECTION  B 

STRAIN  AND  FREQUENCY  DEPENDENT  ATTENUATION  ESTIMATES  IN 
SALT  FROM  SALMON  AND  STERLING  NEAR-FIELD  RECORDINGS 

SUMMARY 

The  amplitude  and  spectra  of  teleseismic  and  regional  seismic  P  waves  depend  on  the 
material  properties  in  the  source  region.  The  effects  of  strain  and  frequency  dependent  attenua¬ 
tion  are  therefore  important  for  understanding  the  source  characteristics  of  coupled  versus 
decoupled  explosions.  Near-field  velocity  and  acceleration  data  from  the  Salmon  and  Sterling 
explosions  in  salt  were  analyzed  for  strain  and  frequency  dependent  attenuation.  Attenuation, 
parameterized  as  1/Q,  was  estimated  for  various  source-receiver  distance  ranges  and  frequency 
bandwidths  by  using  two  methods:  (a)  mean  spectral  ratio  slope  and  (b)  average  amplitude 
spectral  ratio.  Results  from  Salmon  indicate  that  (1)  in  the  frequency  range  of  1-25  Hz,  Q  is 
about  5  to  10;  (2)  Q  increases  with  source-receiver  distance,  suggesting  lower  Q  for  larger 
strain  levels;  (3)  Q  in  the  frequency  range  of  25-50  Hz  is  substantially  higher  than  in  the  fre¬ 
quency  range  of  1-25  Hz,  or  Q  increases  with  frequency.  Analysis  of  data  from  Sterling  indi¬ 
cates  that  Q  is  frequency  independent  and  about  one  order  of  magnitude  larger  than  for  Sal¬ 
mon.  Therefore  the  near-source  attenuation  was  strain  and  frequency  dependent  and  hence 
non-linear  for  Salmon  out  to  a  range  of  600  m. 
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Sterling  at  nearly  equidistant  locations  E6-27  (r  =  318.2  m)  and  E14C- 

36  (r  =  321.0  m). 

5  RVP  spectral  ratios  of  (5A)  Salmon  sensors  E/D,  (5B)  Salmon  E/B,  54 


(5C)  Salmon  sensors  K/H,  and  (5D)  Sterling  sensors  K/H  (pluses)  and 
Salmon  sensors  K/H  (circles).  Note  the  curvature  of  the  Salmon  spec¬ 
tral  ratios  E/D,  E/B  and  K/H  indicating  the  increase  of  Q  with  increas¬ 
ing  frequency.  The  1/Q  values  for  Salmon  K/H  are  indicated  for  the 
1-25,  25-50,  50-100,  and  100-200  Hz  frequency  bands.  Note  the 
difference  in  the  K/H  spectral  ratio  for  common  sensors  pairs  for  Sal¬ 
mon  and  Sterling. 

6  (6A)  1/Q  estimates  for  Salmon  in  the  1-200  Hz  bandwidth  from  the  55 

URH  sensor  pairs.  Mean  UQ  =  0.0244(0.0004),  or  Q  =  40.9(0.7). 

(6B)  1/Q  estimates  for  Sterling  in  the  1-200  Hz  bandwidth  from  the 
URH  sensor  pairs.  Mean  T7Q  =  0.005 1(0.001 1)  or  Q  =  196(40).  (6C) 

1/Q  estimates  for  Salmon  in  the  1-25  Hz  bandwidth  from  the  URH 
sensor  pairs.  Mean  1/Q  =  0.123(0.002),  or  Q  =  8. 1(0.1). 

7  A  weighted  least  squares  fit  of  the  1-25  Hz  Salmon  RVP  spectral  ratio  57 

attenuation  estimates  for  a  fit  against  distance.  1/Q  is  consistent  with  a 

R  decay  with  distance. 

8  1/Q  and  1/Q  estimates  from  RVP  spectral  ratios  as  a  function  of  fre-  59 

quency  bandwidth.  A  smooth  1/Q(f)  curve  is  drawn  through  the  data 
consistent  with  both  the  average  1/Q  values  and  the  slope  of  1/Q  in¬ 
ferred  from  (l/Q-l/Q)/f. 
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INTRODUCTION 

It  is  important  to  understand  the  attenuation  and  other  characteristics  of  rocks  in  the 

' 

vicinity  of  an  explosive  source  because  of  their  direct  impact  on  far-field  amplitudes  and  spec-  i 

i 

tra.  The  detection  of  explosions  and  the  discrimination  of  explosions  from  earthquakes  depends  | 

on  these  far-field  P  wave  spectra,  therefore  in  the  context  of  verification  the  importance  of  J 

understanding  tamped  and  decoupled  explosions  in  salt  becomes  critical  (Hannon,  1985;  ! 

I 

f 

Evernden  et  al.,  1986).  In  this  study,  we  have  analyzed  subsurface  strong-motion  data  from  | 

the  Salmon  explosion  in  salt  and  the  decoupled  shot,  Sterling,  detonated  in  the  cavity  created 
by  Salmon.  Both  of  these  nuclear  explosions  were  recorded  by  a  large  array  of  subsurface 
strong  motion  accelerometers  and  velocity  gauges  (Perret  1968a,  b).  These  data  are  useful  for 
determining  source  differences  between  the  tamped  explosion,  Salmon  and  the  decoupled  shot, 

Sterling,  because  most  of  the  receivers  remained  at  the  same  locations  for  both  explosions. 

The  shot-point  to  receiver  range  covered  by  these  subsurface  instruments  was  from  166  m  to 
740  m,  or  scaled  distances  between  95  m/kt1/3  to  427  m/kt1/3.  The  regional  P  wave  comer  fre¬ 
quency  of  Salmon  suggests  an  equivalent  elastic  radius,  a,  between  240-290  m  (Larson,  1979). 

Data  collected  from  a  variety  of  shots  in  salt  (the  Cowboy  series  of  chemical  explosions,  Sal¬ 
mon,  and  Sterling)  indicate  that  inelastic  effects  extend  to  a  scaled  radius  of  3.6  m/kg1/3  (velo¬ 
cities  of  0.5  m/sec,  strains  of  10-4,  and  a  =  600  m  for  Salmon),  and  that  anelastic  attenuation 
continues  to  at  least  30  m/kg”1/3  (Minster  and  Day,  1986).  Workman  and  Trulio  (1985) 
present  evidence  from  the  Cowboy  Trails  experiments  that  the  inelastic  effects  in  salt  extend  to 
1 1.3  km/kr1/3.  The  decoupled  explosion,  Sterling,  was  shot  in  the  Salmon  cavity  (radius  =  17 
m,  Springer  et  al.,  1968)  and  the  subsurface  measurements  lie  nearly  within  the  linear  region. 

It  would  seem  therefore  that  the  subsurface  measurements  for  Salmon  and  Sterling  were  made 
at  distances  which  may  span  the  transition  from  non-linear  inelastic  into  linear  anelastic 
response  of  the  salt. 
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NEAR-FIELD  DATA  FROM  SALMON  AND  STERLING 

A  detailed  description  of  the  subsurface  and  surface  strong-motion  data  from  Salmon  is 
given  by  Perret  (1968a).  The  azimuthal  relationship  of  each  instrument  boring  to  ground  zero 
(designated  Station  1A)  and  the  range  and  angular  elevation  between  each  instrument  station 
and  the  working  point  are  shown  in  Figure  1  (after  Figure  2.2,  Perret,  1968a). 

The  strong-motion  data  from  Salmon  indicate  peak  radial  particle  velocity  to  attenuate 
with  sour,  e-rcceiver  distance,  r  as  r-1  9  and  the  peak  strains  (peak  velocity/  propagation  velo¬ 
city)  vary  from  about  4xl0~3  to  3xl(T4  (Perret,  1968a,  especially  Figure  5.11).  The  Sterling 
data  indicate  peak  radial  particle  velocity  to  attenuate  approximately  as  r~'  and  the  peak  strains 
span  the  range  of  about  3xl0-5  to  about  7X10-6  (Perret,  1968b,  especially  Figure  3.6).  The 
Salmon  data  have  significant  signal  content  from  1  to  over  200  Flz  whereas  most  of  the  Ster¬ 
ling  data  are  usable  to  400  Hz.  Records  of  all  horizontal  and  vertical  digitized  velocity  gauge 
data  from  Salmon  are  shown  in  Figure  2  (the  12  iastrument  stations  are  designated  as  A,  B,  C, 
D,  E,  F,  G,  H,  J,  K,  L,  and  M).  The  numbers  near  the  observed  waveforms  indicate  zero-to- 
peak  velocity  m/sec.  The  particle  motion  is  up  for  radial  motion  away  from  the  source  or  for 
upward  vertical  motion.  Stations  C,  H,  K,  and  M  have  depths  nearly  equal  to  the  shot  depth  so 
that  vertical  data  from  these  4  stations  represent  transverse  motion.  Similar  available  data  from 
Sterling  are  shown  in  Figure  3. 

Since  the  natural  frequencies  and/or  damping  factors  of  the  recording  instruments  varied 
from  one  location  to  another,  all  available  accelerometer  and  velocity  gauge  data  were  first 
corrected  for  the  instrumental  response  and  velocity  records  were  derived.  Before  comparisons 
between  various  receiver  locations  could  be  made,  it  was  found  necessary  to  correct  the  pani¬ 
cle  velocity  values  for  frequency-dependent  near-field  geometrical  spreading.  Reduced  Velocity 
Potentials  (RVP’s)  were  estimated  from  both  radial  and  vertical  (whenever  not  transverse  to  the 
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Figure  1.  Instrumental  locations  for  the  Salmon  and  Sterling  experiments  (after  Figure  2.2, 
Perret,  1968a). 


Figure  2.  Particle  velocity  records  from  Salmon.  URH  indicates  horizontal  velocity,  UV  verti¬ 
cal  velocity.  Numbers  indicate  peak  velocity  in  m/sec.  Letters  A  through  M  indicate  the 
gauge  locations  used  in  the  analysis. 
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source-receiver  direction)  component  data  from  subsurface  velocity  gauges  and  accelerometers. 
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ESTIMATION  OF  NEAR  FIELD  ATTENUATION 

Although  the  distances  investigated  here  are  admittedly  within  the  non-linear  regime,  we 
have  chosen  to  quantify  the  attenuation  of  the  equivalent  far-field  P-wave  source  (RVP)  with 
the  familiar  parameterization  1/Q.  Therefore,  we  assume  that  at  least  incrementally,  the  spec¬ 
trum  of  the  compressional  pulse  attenuates  as  exp(— =  -^-A(f,r)  is  therefore 

Qc  dr  Qc 

a  quasilinear  approximation  for  the  non-geometric  attenuation  of  the  waveform  spectra.  How¬ 
ever,  in  this  case  we  consider  the  possibility  that  attenuation  is  dependent  on  strain  level  (there¬ 
fore  distance  r)  and  frequency,  Q  =  Q(f,  r).  Work  along  these  lines  has  been  reported  by  Min¬ 
ster  and  Day  (1986)  and  McCartor  and  Wortman  (1985).  Minster  and  Day  (1986)  argued  that 
the  peak  velocity  data  from  the  series  of  experimental  chemical  explosives  detonated  in  salt 
show  that  the  attenuation  must  be  both  frequency  and  strain  dependent.  McCartor  and  Wort- 
man  (1985)  examined  the  Salmon  data  and  concluded  that  the  attenuation  of  the  RVP  must  be 
frequency  dependent.  Previously,  it  was  noted  that  a  frequency  independent  Q  =  10  would 
account  for  the  observed  decay  of  Salmon  peak  velocities  with  distance  like  r-1  9  rather  than  r_1 
expected  from  simple  geometrical  spreading. 

Our  approach  has  been  to  extend  the  analysis  of  1/Q(f)  measured  from  the  Salmon  data 
to  as  many  sensors  available  and  to  investigate  the  strain  dependence,  1/Q(f,  r),  by  examining 
the  distance  ranges  available.  In  order  to  stabilize  the  estimation  of  the  attenuation,  we  have 
estimated  1/Q  for  specific  frequency  bands;  1-25,  25-50,  50-100,  100-200,  and  1-200  Hz. 
These  attenuation  estimates  were  made  in  two  ways;  1)  a  spectral  slope  was  estimated  and  2) 
the  average  spectral  ratio  level  was  estimated  for  each  pair  of  sensors  and  each  frequency  band. 
These  measurements  are  essentially  t*  and  t*  estimates  from  the  spectral  slope  and  spectral 
level  measurements,  respectively.  The  two  estimates  are  related  for  a  frequency  dependent 

—  dt* 

attenuation  by  the  familiar  relation,  t*  =  t*  +  f — —  (  e.g.  Dcr  et  al.,  1982).  For  our  problem  at 

df 
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hand,  the  1/Q  measured  between  two  different  sensors  becomes  the  average  1/Q(r)  between 
two  distance  ranges,  ri(  and  rr  The  attenuation  estimated  from  the  spectral  slope  of  the  spectral 
ratio  between  sensors  i  and  j  is  designated  1/Qjj(f),  while  the  attenuation  derived  from  the  aver- 

_  dO/QTf)) 

age  spectral  level  is  1/Qjj(f).  The  two  are  related  by,  1/Qjj(f)  =  1/Qjj(0  +  f - - - ,  where 


we  assume  that 


1  1 
1/Qij(0  =  — —  [7^—  and  similarly,  T/Q^(f)  =  — — 

rrr i  r,Q(f*r)  rrr*  r,Q(f,r) 


Typical  RVP  spectra  for  Salmon  and  Sterling  are  shown  in  Figure  4.  Note  the  scalloping 
of  the  Sterling  spectra  due  to  the  added  complexities  of  detonation  in  a  cavity,  such  as  pressure 
spikes  (Larson,  1979).  Typical  RVP  spectral  ratios  are  shown  in  Figures  5 A,  5B,  5C  for  Sal¬ 
mon.  Note  the  curvature  of  the  plotted  spectral  ratios.  It  is  not  possible  to  produce  such  a 
spectral  ratio  without  an  increase  in  Q  with  increasing  frequency.  This  is  the  same  basic  meas¬ 
urement  used  by  McCartor  and  Wortman  to  estimate  frequency  dependent  Q(f)  for  the  Salmon 
data.  If  we  compare  the  spectral  ratios  from  Salmon  to  those  of  Sterling  in  Figures  5C  and 
5D,  we  see  that  although  the  Sterling  data  is  modulated  by  spectral  scalloping,  the  average  Q 
derived  from  Salmon  cannot  accommodate  the  data  from  Sterling.  The  average  Q  from  the 
Sterling  data  in  the  1  to  200  Hz  bandwidth  is  200  while  the  average  Salmon  Q  derived  from 
both  spectral  slopes  and  absolute  spectral  levels  is  less  than  40.  These  average  1/Q  estimates 
are  distinct  at  the  three  a  confidence  level  (see  Table  1).  Previous  estimates  of  the  near-field 
Sterling  attenuation  were  Qs  =  35  and  QP  =  70  by  Langston  (1983)  who  modeled  the  rise 
times  of  P  and  S  acceleration  waveforms  as  a  function  of  distance  in  the  salt. 


The  1/Qjj  and  1/Qjj  estimates  and  their  estimated  uncertainties  were  tabulated  for  the  1-25, 
25-50,  50-100,  100-200,  and  1-200  Hz  bandwidths.  The  estimates  for  the  1-25  and  1-200  Hz 
bands  can  be  seen  in  Figures  6A,  6B,  and  6C  where  the  horizontal  bars  represent  the  distance 
range  over  which  the  attenuation  estimates  were  made.  Weighted  averages  over  all  sensor 
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Figure  4.  Spectra  of  RVP  obtained  from  particle  velocity  data  from  Salmon  and  Sterling  at 
nearly  equidistant  locations  E6-27  (r  =  318.2  m)  and  E14C-36  (r  =  321.0  m) 


Figure  5.  RVP  spectral  ratios  of  (5A)  Salmon  sensors  E/D,  (5B)  Salmon  E/B,  (5C)  Salmon 
sensors  K/H,  and  (5D)  Sterling  sensors  K/H  (pluses)  and  Salmon  sensors  K/H  (circles).  Note 
the  curvature  of  the  Salmon  spectral  ratios  E/D,  E/B  and  K/H  indicating  the  increase  of  Q  with 
increasing  frequency.  The  1/Q  values  for  Salmon  K/H  are  indicated  for  the  1-25,  25-50,  50- 
100,  and  100-200  Hz  frequency  bands.  Note  the  difference  in  the  K/H  spectral  ratio  for  com¬ 
mon  sensors  pairs  for  Salmon  and  Sterling. 


SALMON  UR  1.0  -  200.0  H* 


1/(2  -  0.0244(0.0004) 


00  26.00  36.00  ‘♦6.00  56.00 

DISTANCE  *10' 


66.00  76.00 


STERLING  UR  1.0  -  200.0  Hi 


<k  1/(2  -  0.0051(0.0011) 


,00  26.00  36.00  46.00  56.00  66.00  76.00 

DISTANCE  *10' 

SALMON  UR  1.0  -  25.0  Hi  C 


ik  1/(2  -  0.123(0.002) 


16.00  26.00  36.00  46.00  56.00  66.  00  76.  GO 

DISTANCE  *10' 

Figure  6.  (6A)  1/Q  estimates  for  Salmon  in  the  1-200  Hz  bandwidth  from  the  URH  sensor 
pairs.  Mean  1/Q  =  0.0244(0.0004),  or  Q  =  40.9(0.7).  (6B)  1/Q  estimates  for  Sterling  in  the 
1-200  Hz  bandwidth  from  the  URH  sensor  pairs.  Mean  1/Q  =  0.0051(0.001 1)  or  Q  =  19b(40) 
(6C)  1/Q  estimates  for  Salmon  in  the  1-25  Hz  bandwidth  from  the  URH  sensor  pairs.  Mean 
T7Q  =  0.123(0.002),  or  Q  =  8. 1(0.1). 
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combinations  are  presented  in  Tables  2  and  3.  The  average  attenuation  estimates  from  the  hor¬ 
izontal  velocity  (URH),  vertical  velocity  (UV),  horizontal  acceleration  (ARH),  and  vertical 
acceleration  (AV)  sensors  are  listed  separately  as  well  as  all  sensors  together  (ALL)  in  Tables 


2  and  3.  Q(f)  must  substantially  increase  with  increasing  frequency  for  Salmon  between  the 
1-25  and  100-200  Hz  bands. 


TABLE  1.  COMPARISON  OF  1/Q  FOR  SALMON  AND  STERLING 


SALMON 

STERLING 

FREQ  BAND 

1/Q 

a 

Q 

c 

n 

1/Q 

a 

Q 

a 

n 

1-25 

0.1235 

0.0014 

8.1 

0.1 

99 

0.0026 

0.0002 

389 

37 

87 

1-200 

0.0290 

0.0002 

34.5 

0.3 

99 

0.0026 

0.0005 

384 

71 

87 

TABLE  2.  MEAN  1/Q  ESTIMATES  FOR  SALMON 


FREQ  BAND  ALL 


1-25 

25-50 

50-100 

100-200 

1-200 


0.16(0.08) 

0.096(0.006) 

0.059(0.004) 

0.049(0.003) 

0.049(0.003) 


AV 

0.21(0.02) 

013(0.01) 


ARH 

0.09(0.04)  0.14(0.02) 

0.09(0.01)  0.09(0.01) 

0.056(0.004)  0.055(0.004) 

0.0042(0.003)  0.037(0.003) 

0.0047(0.005)  0.044(0.003) 


0.073(0.007)  0.056(0.004) 

0.051(0.003)  0.0042(0.003) 

0.061(0.004)  0.0047(0.005) 


URH 

0.15(0.01) 

0.10(0.01) 

0.059(0.007) 

0.040(0.006) 

0.0047(0.007) 


TABLE  3.  MEAN  1/Q  ESTIMATES  FOR  SALMON 


FREQ  BAND  ALL 


1-25 

25-50 

50-100 

100-200 

1-200 


0.09(0.01) 

0.009(0.008) 

0.026(0.002) 

0.026(0.002) 

0.029(0.002) 


AV 

0.18(0.04) 

0.07(0.02) 


_ ARH 

0.04(0.02) 

-0.02(0.01) 


0.019(0.003)  0.026(0.002) 

0.034(0.005)  0.026(0.002) 

0.033(0.003)  0.032(0.002) 


UV 

0.11(0.02) 

0.01(0.02) 

0.15(0.007) 

0.025(0.0004) 

0.025(0.004) 


URH 

0.119(0.006) 

0.01(0.01) 

0.026(0.008) 

0.016(0.005) 

0.026(0.006) 


The  1/Q,j  data  were  inverted  for  distance  dependence,  1/Q(r),  by  a  weighted  least  squares 
procedure.  The  average  1/Q  was  determined  for  each  distance  range  bracketed  by  the  sensor 
distances.  The  results  are  shown  in  Figure  7  for  1-25  Hz.  The  Salmon  1-25  Hz  near-held 
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attenuation  declines  by  nearly  50%  between  200  and  600  meters.  Because  the  attenuation  is 
smaller  for  the  higher  frequencies,  it  is  less  certain  that  there  exists  a  distance  dependence  in 
the  rate  of  attenuation  for  the  25-200  Hz  range. 

The  1  to  25  Hz  1/Q(r)  appears  to  decline  like  r-1/2  rather  than  like  r'1 9  as  would  be 
expected  if  the  attenuation  rate  were  proportional  to  the  strain  level,  £.  The  near-field  attenua¬ 
tion  of  the  Sterling  RVP  was  so  little  that  no  reliable  distance  dependence  can  be  seen  over 
such  a  short  distance.  Therefore  it  may  be  possible  to  coasider  the  Sterling  Q  as  the  low 
strain-rate  background  Q  level  for  salt.  It  is  interesting  that  McCartor  and  Wortman  (1985,  see 
appendix  A)  have  shown  that  the  peak  velocity  data  for  explosions  in  salt  could  be  consistent 
with  attenuation  proportional  to  Ve  or  to  £2,  as  well  as  directly  proportional  to  £.  Since  the 
Minster  and  Day  model  was  derived  from  data  at  strains  less  than  the  Salmon  data  it  is  not 
clear  whether  a  transition  from  £1/2  to  £  occurs  near  strains  of  10“ 4  or  whether  another  parame¬ 
terization  of  the  Cowboy  data  is  possible. 

The  data  require  that  the  attenuation  rate  be  frequency  dependent  as  well  as  depend  upon 
the  strain  level.  The  net  difference  between  the  average  1/Q  and  the  1/Q  in  the  1-25  Hz  range 
is  consistent  with  an  average  slope  of  1/Q(f)  of  -0.002  Hz-1  in  the  1-25  Hz  bandwidth.  This 
is  also  consistent  with  the  net  average  difference  between  the  average  1/Q  and  1/Q  between  1- 
25  and  25-50  Hz  bandwidths;  -0.002  Hz-1.  These  data  are  summarized  in  Table  4  for  each 
frequency  band.  Assuming  that  1/Q,  1/Q  and  (l/Q-l/Q)/f  are  smooth  functions  of  frequency, 
we  can  derive  a  model  that  accommodates  both  the  absolute  spectral  ratio  levels  as  well  as  the 
average  slopes  of  the  spectral  ratios.  This  1/Q(f)  model  is  plotted  in  Figure  8.  The  model 
differs  from  McCartor  and  Wortman’s  model  where  they  parameterize  their  model  as  1/Q(0  = 

Yo  +  Yi  r1.  fAiffl  =  a  constant  for  the  McCartor  and  Wortman  parametrization.  We  find 
df 

that  this  parameterization  is  inconsistent  with  the  data  of  Tables  2  and  3.  The  changing  curva- 
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NEAR-FIELD  SALMON  ATTENUATION  ESTIMATES 
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Figure  8. 
bandwidth. 


1/Q  and  1/Q  estimates  from  RVP  spectral  ratios  as  a  function  of  frequency 
A  smooth  1/Q(f)  curve  is  drawn  through  the  data  consistent  with  both  the  average 


1/Q  values  and  the  slope  of  1/Q  inferred  from  ( 1/Q —  l/Q)/f. 
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DISCUSSION  AND  CONCLUSIONS 

We  see  that  the  Salmon  near-field  data  indicates  that  the  non-linear  near-field  Q(f,  r) 
increased  with  increasing  frequency  and  decreasing  strain  rate.  The  Sterling  linear  near-field  Q 
was  nearly  frequency  independent  and  greater  by  an  order  of  magnitude.  However,  to  clearly 
see  these  relationships  it  is  necessary  to  average  over  either  distance  ranges  or  frequency  bands 
to  attain  the  respective  resolution  of  1/Q(f,  r)  in  either  the  the  frequency  or  spatial  domains. 
While  the  Sterling  data  are  consistent  with  the  theory  of  Minster  and  Day  (1986)  based  on  the 
Cowboy  chemical  explosions  in  Salt,  the  Salmon  1/Q  declines  at  a  slower  rate,  =1/Vr,  than 
would  be  expected  if  the  attenuation  rate  were  directly  proportional  to  strain.  These  two 
nuclear  explosions  bracket  the  Cowboy  chemical  explosion  data  modeled  by  Minster  and  Day. 
Another  possibility  may  be  consistent  with  the  results  of  Figure  7,  that  the  attenuation  rate 
changes  near  300  meters  and  is  nearly  constant  for  ranges  greater  than  300  meters.  This  is 
roughly  coincident  with  the  equivalent  elastic  radius  inferred  from  the  regional  P-wave  comer 
frequency. 

We  see  that  explosions  decoupled  in  salt  may  be  enhanced  in  high  frequencies  over  a 
well  coupled  shot  in  salt  because  of  the  strain  and  frequency  dependence  of  near-field  attenua¬ 
tion.  Springer  et  al  (1968)  estimated  the  decoupling  factor  of  Sterling  to  be  70(20)  at  1  to  2 
Hz  using  surface  seismic  recording  stations.  Blandford  and  Woolson  (1979)  used  regional  data 
to  estimate  the  Salmon/Sterling  spectral  ratio  for  P  waves  to  be  about  1000  at  1  Hz  and  only 
about  17  at  25  Hz.  The  decline  of  decoupling  at  high  frequencies  becomes  a  critical  factor  for 
detection  and  discrimination  of  decoupled  shots.  In  case  of  the  Salmon/Sterling  experiments, 
the  differences  in  near-field  attenuation  may  account  for  as  much  as  a  factor  of  2  at  10  Hz,  and 
5  at  100  Hz,  if  the  proposed  Q(f,r)  model  is  valid  in  the  range  from  200  to  1000  meters. 
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